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Abstract. We study the internal circulation within the cocoon carved out 
by a relativistic jet emanating from an AGN, first developing model and then 
validating it using a series of numerical simulations. We notice that a significant 
increase of circulation in this flow arises because gradients in the density and en- 
tropy develop near the hot spot, as a consequence of Crocco's vorticity theorem. 
We find simple solutions for the streamlines, and we use them to predict the 
mass inflow rates towards the central regions. The 2D simulations we perform 
span a rather wide range of mechanical jet's input power and Black Hole masses, 
and we show that the predicted nuclear mass inflows are in good agreement with 
the theoretical model. 

We thus suggest that these backflows could (at least partially) feed the AGN, 
and provide a self-regulatory mechanism of AGN activity, that is not directly 
controlled by, but possibly controls, the star formation rate within the central 
circumnuclear disk. 



1. Introduction 

The presence of compact objects (hereafter Black Holes, BHs) in the central 
dense regions og galaxies is often related to the AGN phenomenon . More- 
over, since the original theoretical suggestion by ISilk &: Reed (jl998l ). signifi- 
cant observational evidence has been accumulated concerning the impact of nu- 
clear activity on the global stellar evolution within the host galaxy (see e.g. 
Schawinski et alJl2006lpYi et al.1 [20071 : iKavirai et al.ll2007l . [20081 . for some recent 
work). In addition to its effect on global star formation, the presence of an AGN 

seems also to be conn ected with circumnuclear starbursts on small (from -parsec 

to kilq-parsec) scales ( Storchi-Bergmann et al.ll200ll ; iGonzalez Delgado Sz Cid Fernandes 



20051 : ISarzi et al.l [20071 : iDavies et al.ll2007l ~ However, despite all this observa- 



tional evidence, the connection between AGN activity, negative/positive feed- 
back on star formation, and local starbursts is not well understood. 
If the main physical agent for this connection is the interaction between the 
jet and the host galaxy's interstellar medium (hereafter ISM), then it is impor- 
tant to model the physical consequences of the propagation of AGN relativistic 
jets into the ISM, and the mechanism which feeds the parsec-scale accretion 
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disc around the central BHs. R ecent simulations ([Sutherland fc Bicknelll 120071 ; 
Antonuccio-Delogu fe Silkll2008l ) have begun to self-consistently model the feed- 



back of a relativistic jet on its host galaxy, a nd the global conseq uences for e.g. 
blue-to-red cloud migration and downsizing ( Tortora et al J [20091 ). Here we will 
focus our attention on the large scale structure of the flow within the cocoon, 
and on its consequences for the accretion onto the central BH. 



2. Formation of the backflow 

We distinguish 3 different sections in the circulation (see Fig. [1]): propagation 
near the jet, then from the hotspot along the inner part of the bow shock, and 
along the meridional plane. 



2.1. Model 



Falld (|!99lh and H aiser Alexander] (jl997l ) have shown that a recollimation 
shock (hereafter RS) forms at some distance along the path of the jet. The 
post-shocked gas then accumulates into a region confined between the RS and 
the outer surface of the bow shock (see Figure [1]), with a density n^g ~ 7nj, as 
appropriate to shocked, hot (relativistic) gas. We generically call this region the 
hot spot (HS). 

Most of the cocoon is occupied by very low density, high temperature gas, 




Figure 1. Schematic model of the origin of circulation near the hot spot 
and in the equatorial plane of the cocoon. We assume that the hotspot (HS) 
is bounded by a lateral spherical section having curvature radius r s . 



which keeps the jet pressure confined and collimated, even in the absence of 
magnetic fields. In this region the gas is in a turbulent state, with very little or 
no systematic motions. This region however does not extend directly down to 
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the boundary of the jet: the initial shear of the latter induces parallel motions 
of the gas near and outside the jet axis. This gas eventually reaches the HS, 
post-shocked region, where the density is significantly higher than in the cocoon. 
Obviously, the shocked gas within the hot spot has a higher entropy than that in 
the cocoon. Under these conditions, Crocco's theorem states that a circulation 



arises within a compressible fluid, even if in laminar motion (see ICapl 12001 
for a more recent discussion). In planar geometry the only component of the 
circulation is normal to the flow plane, and its magnitude is given by: 

ds dh 

uv = T— - — (1) 

an an 

where: n is the normal to the direction perpendicular to the streamline, co = | ui \ 
is the modulus of the circulation, and s and h$ = h + v 2 /2 are the specific en- 
tropy and stagnation enthalpy, respectively. 



2.2. Flow near the Hot Spot 

We will approximate the recollimation shock surface before the hotspot as pla- 
nar, thus the entropy is constant across the streamline, and the flow will not 
gain any circulation. However, the gas flowing near and outside the jet, will 
also eventually reach the HS: the boundary layer separating this region from 
the turbulent region will be a curved surface, joining the recollimation shock 
to the bow shock (see Fig. [TJ. We then assume that this off-axis flow motion 
is predominantly directed along the z— direction, thus: v = v u \ z i. Hereafter, 
subscripts u and d will denote quantities computed in the upstream and down- 
stream regions, respectively, where by the latter we mean the HS, as shown in 
Figure [TJ In order to compute Uhs, the circulation in the hot spot, we will make 
use of the expression for the circulation near a curved shock, derived bv lShapirol 
(|!958l . eq. 8): 

2v z \ u {Ml -I) | cos 0| 
Uhs r,( 7 + l)M2[2 + ( 7 -l)M2] U 

where: M n = (v z /c s ) cos# is the Mach number in the upstream flow region. As 
we show in Antonuccio-Delogu and Silk (2009, submitted) one can obtain an 
exact solution in 2D spherical coordinates for the velocity field after the shocked 
layer: 

( 4v z \ u {Ml -I) | cos 0| 

^ j "(7 + l)M2[2 + ( 7 -l)M2] [6) 

where <p is a polar angle and. We see that the azimuthal component of velocity 
after the shock becomes negative, and the polar component in the downstream 
region acquires a positive value: thus, a backflow develops. 

2.3. The flow along the bow shock 

After the HS, the backflow enters into the inner part of the bow shock (Figured]), 
and the magnitude of the rotation is determined by the conservation law for 
vorticity in 2D: 

I u I . . 
= const (4) 
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Figure 2. Left: The excursion of | v d \ z /v u \ z |, as measured in eq. G2 for 
different values of the upstream Mach number. The flows within the cocoon 
usually are transonic, with Moo —1 — 2. Right:Thc scaled density within 
the bow shock as a function of a/ Act, the semi- major axis measured in units 
of the width of the spheroidal shell containing the shock. We plot the ratio 
p/(3poflo/(^ a ) 2 )' f° r t wo different values of the aspect ratio R. 

Thus, in order to determine | co |, we have to determine the density inside the bow 
shock. Assuming that all the matter which was inside the cocoon is compressed 
within the bow shock, the final result reads: 



where po is the central density of the galaxy halo, and we have defined: A = 
a^/Aa, where aj,Aa are the semimajor axis and the width of the bow shock, 
respectively. Note that the time-dependance of the density enters through the 
semi-major axis (a = a(t)). Integrating eq. [5] we obtain the velocity along the 
bow shock, vf. the result is shown in Fig. [3] {left). We again notice that the 
maximum variations are always around unity, along the streamlines, and the 
velocity tends to decrease or to stay almost constant, as we approach the polar 
regions. 

In summary, our model makes a few predictions. First, at least during the initial 
phases of the expansion, a circulation arises inside the cocoon, induced by the 
presence of a high density and entropy region (the Hot Spot). The backflow 
which develops tends to follow the bow shock, and then bends back near the 
meridional axis, perpendicular to the jet. For symmetry reasons, this flow will 
converge back towards the jet. 

2.4. Numerical simulations 

The numerical simulations confirm the numerical model sketched (Fig. [31 right). 
It is interesting to observe that the agreement with the model predictions is 
quite satisfactory, for a wide range of central densities and jet's injection powers. 
Note that the in these simulations the main parameters like central halo mass, 
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Figure 3. (Left): Velocity along the bow shock, for different aspect ra- 
tios. (Right): Mass inflow rates around a central disk region. The continous 
lines are for simulations s(100,200,300)av, dotted lines for s(100,200,300)pl, 
dashed lines for s(100,200,300)ml. The dashed lines are the fits from the 
model given in the text. The red vertical lines mark the epoch when the 
reconfincmcnt shock is destroyed and the jet starts to propagate freely, being 
only confined by the cocoon's ram pressure. 




Figure 4. (Left): Velocity field before the destruction of the jet's shockby 
KH instabilities. The general pattern is well described by the model sketched 
in this paper. (Right): pattern of the flow after the disappearence of the shock 
along the jet. Note that the colour scale (not shown) is much smaller than in 
the figure on the right, thus density differences are in reality very small. 



BH mass, jet's injection power, are chosen according to scaling relations (see 
Antonuccio-Delogu and Silk, 2008, for details). 

The simulations also show that two different circulation regimes describe the 
flow within the cocoon (see Fig. [3]). The presence of the shock before the HS 
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determines a flow similar to that described in our model. However, after a 
while this shock is destroyed by Kelvin-Helmholtz instabilities, and the hotspot 
propagates freely (Fig. UJ right). The gas near the jet is now reflected by the 
HS, and determines a counterstreaming flow which propagates back towards the 
central disc. The net mass inflow rate is however much less. 

3. Conclusion 

Realistic simulations of the flow pattern within an expanding cocoon show that 
a backflow develops soon after the cocoon forms. This determines a coherent 
flow towards the central parts of the AGN, whose magnitude is of the order 
of few Mq yr _1 , i.e. the right order to feed the AGN and sustain its activity. 
A more detailed investigation of the consequences of this backflow is presented 
elsewhere (Antonuccio-Delogu and Silk, 2009, submitted). 
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